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SHORT OPTICAL PULSES EVOLUTION 
IN NONLINEAR DIELECTRIC THIN LAYERS 


Mihai C. PISCUREANU!, Vicentiu CONSTANTIN? 


Rezumat. fn lucrare se analizeaza propagarea impulsurilor temporale intr-un strat 
dielectric subfire neliniar constand dintr-un material dispersiv neomogen, in pre- 
zenta efectelor de rezonanta. Am aplicat o metoda perturbationala si cu ajutorul 
acesteia am obtinut solutii solitonice de tip anvelopa, stabilind astfel conditiile de 
existenfa a solitonilor strdlucitori si respectiv intunecafi cu ajutorul coeficientilor 
ecuatiei Schrédinger neliniare si al relatiei de dispersie a undei cuasiplane. Aceste 
rezultate teoretice au fost aplicate in cazul unei structuri dielectrice planare cu 
SiO». In acest caz am ardtat ca in vecinatatea frecventei de rezonanfa si a frecventei 
de tdiere superioare solitonii inceteaza de a mai fi impulsuri scurte temporal. 


Abstract. In this paper we have analyzed the propagation of short temporal pulses 
in a nonlinear dielectric thin layer consisting of an inhomogeneous dispersive 
material in the presence of the resonance effects. We used a perturbational method 
and we obtained envelope solitons. Thus, we have established the existence 
conditions for the bright and dark solitons from the nonlinear Schrodinger equation 
coefficients by using a cuasiplanar wave dispersion relation. The theoretical results 
have been applied in the case of a typical SiO, slab. Consequently, we have shown 
that in the vicinity of the resonance frequency and the upper cut-off frequency, the 
solitons cease to be temporally sharp pulses. 
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1. Introduction 


In this paper we analyze the propagation of the pulses in a dispersive Kerr-type 
medium, represented by an infinitely long dielectric, without losses and with a 
parabolic profile of the refractive index. Our analysis has the following stages: 


a. The effects of the transversal inhomogeneities are decoupled from the 
longitudinal propagation in the planar dielectric medium. This is achieved by an 
averaging method, thus obtaining in the monomodal case a nonlinear wave 
equation describing the longitudinal propagation of the pulses. 


b. We consider all the harmonics generated by the nonlinearity by 
expanding the unknown function in a power series of the ¢ parameter (¢<<1). 
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We analyze the case of the temporally short pulses, for which the parameter ¢ is 
given by the ratio between the spectral width Aq@ and the carrier angular 
frequency @. 


c. We establish the integrability conditions for the order of ¢, thus 
obtaining the corresponding group velocity expression and the dispersion 
relations. 


From the integrability conditions corresponding to o(.3) we establish the 
nonlinear Schrédinger equation (NLS). 


By means of NLS we establish the conditions leading to the generation of bright 
and dark solitons in the presence of the medium dispersion. 


2. Theoretical model 


a. The refractive index n(x), where x is the transverse coordinate, is 
given by: 


x =2 
2 1-24{ 2) + K.|E| (1) 


where mj is the linear refractive index for x=0, A=(n, —109)/ny, is the relative 


variation of the film refractive index with respect to that of the coating; K, is the 


e 
Kerr coefficient and E is the electric field intensity vector. 


Expression (1) is valid as long as the resonant effects can be neglected. 
The resonance effects denoted by AP;, can be described as a linear contribution to 


the polarization vector Py, =€o 122() ~ifE, and are given by a convolution 
integral: 


aE 
AP, =é 9 - fa’ [ar o(z — z')- (z= 23t-2)E(&, z',t') (2) 
—0 —00 
&9 is the dielectric constant of the vacuum, 0 (e) is the delta function, expressing 
the spatial homogeneity along z-axis, and b(z,t) is the anomalous part of the 
response function of the dispersive medium. 
Considering the even TE mode case with E = jE y (x Zz, t) and using the Maxwell 


equations, we get the following wave equation: 
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(3) 


where AP; = |AP_| : 


In the linear case and in the absence of the resonance effects (K e =0,AP; = 0), 
the solution of equation (3) can be expressed as Hermite polynomials depending 
of the transverse coordinate x. For the linear monomode case, E y becomes: 


Ey(x,z,t)=Ep exp(- x? 1 2), exp|i(@ot — Fz) (4) 


where Ep is the normalizing constant, @g is the angular frequency, w is the 


beam radius and fp is the wave vector in the considered medium. 
In the nonlinear case, the effects of the transverse inhomogeneities are taken into 
account by averaging. Thus, taking E,, (x, Z, t) as 
Ey(x%,zt)=U(z 1)-exp{- x2 /w?} (5) 
y ue ? ) 


after applying an averaging method, we get the following onedimensional 
nonlinear wave equation: 


eu um f=) K, “(op | ‘). 
Oe- w* A SAN) Be? Wa 3c? ar? 

\ Oe? ; 
NS are 


(6) 
[ae fades KG z'). D(z (4 F t')-U(z08) 

—00 —00 

b. The perturbational method 


Suppose that the unknown function U (z,t) can be expanded by the powers of the 
€ parameter (e << 1): 


U(z,t) = >, > e"UF (E,0)-exp[-il(kz — at] (7) 


l=-oa=1 
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where k is the wave number, @ is the angular frequency and the summation is 
made for all the harmonics generated by the nonlinear response of the 


polarization. U : is the envelope of the /-th harmonics of the approximation of 


order a, slowly varying with z and tf. 


The variables € and 7 are defined by the following scaling relations: 


E= Ez, | (8) 
Veg 


where v, is the group velocity. 


From relations (6) — (8) we get the following equation for the /-th component 


U; (é,2): 


4 0 ¢ A ne é Z 
a 2 . a a 1 . a 
> E s° ———-—ilk| U RK —t.= 1- s—t+ilom | Ur - 
! w- ae o)( Ot 5 / 


2 
K. = = b+c 0 a b c 
Wee De 2 . [eZ ito (us fe Wer | 
m=—o c=1 
1 = (-i) “ph 0 ) G\ Us De ep 
_ = —S_— — sii eet 
p> aim? (Sen!) | araem || Baran” Jc 
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(9) 
where p=(VA/ 2\(kond) and by, is the Laplace transform of b(z,1) 
function). 


In order to separate different orders of ¢, we introduce the following operator: 


Qo 1 
L=-k? +——|n? (1 Begs (10) 
Cc WwW 


where by =bo,x- 
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Finally we get: 
- third order: 
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m=—oo 
(11) 
c. Integrability conditions 
From (11) we obtain the integrability conditions for O (<3 ) 
ve) (12) 
on {1 k)é cae 
+h ibe nh \hS Pheoppr| = 0 13 
do 2 al a Uap) Pe ,” 
1a°L otk OU; fal, w 
3 =% Ni r n2(1- p) +b | i 
260" dw} dr* | OK\ ce? 
(14) 
2 au! K 
BN Ft SU WU pom =O 
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Equation (14) is a nonlinear Schrédinger-like equation, whose analysis is useful in 
establishing the conditions for the generation of bright and dark solitons. 


By means of relations (10) and (14) for one mode propagation, we get: 
2 eu! 24 6°U} 3K, 
(2 2 1, 2k OU; 3Ko 


—k 
2 Ok 0g aa” ar? Oe 2 


ut}  =0 (15) 
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When the resonance effects can be neglected, (be = 0) , equation (15) becomes: 


1 2 2771 
1 OU 2 
ra ae + wit =0 (16) 
OC 2G" 
where: 
3K 
Y= < 7 (17) 
2./2¢ 
The relation (15) can be written under the form: 
OK 
K—_~ 2 D 2 
FY wet 2 by CP] Reme, Bo 
2 eK 2 OK 
where 
>=O0,, k=". Ket SOU, aX, cr=T-0;' > 
G a) 


@, is the resonance frequency of the function b,,, and Ug is the initial amplitude 
of the pulse. 


From relation (18) we can establish the conditions for the generation of bright and 
dark solitons, respectively. 


In order to characterize the pulse width we introduce the following parameter: 
=1/2 
2, 2, 
at : si k ake | (20) 
29) 


named characteristic frequency, and being correlated to the temporal width of the 
pulse. 


In the absence of the resonance effects (be = 0) because v is always positive and 


2 
lek. ti [2] bis 


2 da’ (20,0 O, 


is negative (v, 20; >0) , it results that the equation (16) has only bright solitons 


solutions. 
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3. Analysis of the bright and dark solitons generation conditions in a planar 
dielectrical structure 
Suppose that the function b(z,t) is given by: 


2 
Op. 
b(z,t) = b(t) = ——sina,t (22) 


ée 


and consequently the function b,, , is of the form: 
big. = By => (23) 


In this case, b,, corresponds to the real part of the electric susceptibility and a, 
and @, are the electric dipole resonance frequency and the electronic plasma 


frequency, respectively. 


By means of the expression of b,,, we can obtain the dispersion relation L=0 
under the form: 


@O 
7 oe ol P (24) 


In order to calculate the characteristic frequency © and establish its sign we use 
the relation (24) under the form: 


=e la? - 0? 0? Mo? | 
n; (L— p) \ (9?) ay 


where Q, and Q; are the upper and lower cut-off frequencies respectively: 


2 
oO 
20¢n; (=?) 5g 
: (26) 
5 Cc 
OT a a5 


VOew Ny (l- p) 


For the cuasiplanar wave approximation, the propagation is possible when 
OFe Os 1-or £2582. 
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When resonance is neglected, the waves propagate when Q; <Q, the dispersion 
relation becoming: 


/2 
K(Q)=+4n,(1— p)/2(a2 - 9?) (27) 
and due to the fact the upper cut-off frequency becomes equal to 1, we get: 
I Q? 
= nf =p} — (28) 
OQ ORG; 


Because Q, is always negative it follows that equation (18) will have only bright 
soliton-like solutions. 


In figures 1 and 2 the characteristic frequency ©, is represented as a function of 
the normalized frequency © and of the normalized wave vector K respectively. 


(10,1) 
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Fig. 1. Characteristic frequency Q dependence on normalized frequency Q for bright and dark 
solitons case. 

The positive values (dark solitons) and the negative values (bright solitons) of Qo 
are shown by the fine dashed lines and by solid lines respectively. 

The long dashed lines refer to the case of the resonanceless dispersion; in this case 
the absolute value of © behaves like a linear function of © and of K 


respectively for values of © much larger than the cut-off frequency 
Q; (Q, =1). 
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If the resonance effects are present, then we must take into consideration the 
upper cut-off frequency Q, (Q ae 1) and the resonance frequency Q=1. 


The graphs have been drawn for planar SiOz media, the parameters being properly 
chosen to permit the propagation of a single mode, namely: dielectric film 
thickness d = 2 um; the relative variation of the core refractive index with respect 


to that of the cover, A = 0,01; the vacuum wavelength, 27/ky = 0,588 yam. 


The electric dipole resonance frequency is 20358 THz and the electron plasma 
frequency is 21321 THz. 


The corresponding normalized cut-off frequencies are ©); =8,62-10°° and 
Q,, = 1,23 respectively. 
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Fig. 2. Characteristic frequency 9 dependence on normalized wave vector K for bright and dark 
solitons case. 


From figures 1 and 2 we noticed the existence of a normalized transition 
frequency ©,, around which the characteristic parameter 9 has high values, 
which implies the fact that solitons with short temporal width propagate in the 
vicinity of that frequency. 


It is significant the fact that nature of solitons changes when passing through the 
transition frequency, namely from bright solitons to dark ones. 


Dark solitons propagation takes place within the frequency range between the 
transition frequency Q, and the resonance one Q=1. 


68 Mihai C. Piscureanu, Vicentiu Constantin 


The temporal width of the temporal “kink’-like solitons increases far from the 
transition frequency ©, and rapidly decreases when © approaches the resonance 
value Q=1. 


On the other side, the bright solitons have narrow width only for frequencies 
much larger than the cut-off frequencies Q; and Q,. 


Conclusions 


It has been shown_that in the absence of the resonance effects in the planar 
waveguide material only bright soliton-like solutions exist. 


For short temporal pulses a bright-to-dark solitons transition occurs at a critical 
frequency Q,. 


However, at the frequency Q, (Q, <Q) or at resonance (Q <1) the bright and 


dark soliton solutions are not short temporal pulses any longer, due to the rapid 
decrease of their characteristic frequency Qo. 


In conclusion, the existence of a bright-to-dark solitons transition at the frequency 
Q,, lying between the lower cut-off frequency ©; and the upper one Q, is a 
characteristic feature of the short temporal pulses evolution and it is due to the 
presence of the resonant dispersive properties of the dielectric medium. 
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